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Abstract 

o: 

O. 
O ■ 

p^ . Interferometric gravitational wave detectors operate by sensing the differ- 

Oj. ential light travel time between free test masses. Correspondingly, they are 

^ \ sensitive to anything that changes the physical distance between the test 

masses, including physical motion of the masses themselves. In ground-based 

detectors the test masses are suspended as pendula, in order that they be 

approximately "free" above the pendulumn frequency. Still, thermal or other 

excitations of the suspension wires' violin modes do impart a force on the 

masses that appears as a strong, albeit narrow-band, "signal" in the detec- 

^p ' tors wave-band. Gravitational waves, on the other hand, change the distance 

(^ ■ between the test masses without disturbing the suspensions. Consequently, 

^^ ■ violin modes can confound attempts to observe gravitational waves since "sig- 

-■««« , nals" that are correlated with a disturbance of the suspension violin modes 

qh| are not likely due to a passing gravitational wave. 

^ \ Here we describe the design of a Kalman filter that determines the time- 

p|J' dependent vibrational state of a detector's suspension "violin" modes from 

^ ■ time dependent observations of the detector output. From the estimated 

^ . state we can predict that component of the detector output due to suspen- 

H ' sion excitations, thermal or otherwise. The wire state can be examined for 

■ ■ ' evidence of suspension disturbances that might, in the absence of such a di- 

agnostic, be mistaken for gravitational wave signals. Additionally, from the 
wire state we can subtractively remove the contribution from suspension dis- 
turbances, thermal or otherwise, from the detector output, leaving a residual 
free from this instrumental artifact. We demonstrate the filter's effectiveness 
both through numerical simulations and application to real data taken on the 
LIGO 40 M prototype detector. 
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I. INTRODUCTION 

Within the next several years five kilometer-scale interferometric gravitational wave de- 
tectors (GEO [1], LIGO [2] and VIRGO [3]) will begin searching for astronomical sources of 
gravitational waves, something the smaller TAMA [4] prototype detector has already begun. 
The fundamental limitation on the sensitivity of these terrestrial detectors is seismic noise 
at low frequencies, thermal excitations of the mirror substrates and their suspensions at 
intermediate frequencies, and laser shot noise at higher frequencies. Here we focus on that 
part of the suspension noise associated with the violin modes of the suspension wires. 

Gravitational waves interact with interferometric detectors by changing the separation, 
and thus the round-trip light travel time, between the mirrors in the detector arms. In order 
to be sensitive to passing gravitational radiation, the mirrors in terrestrial interferometeric 
detectors are suspended as pendula in order that they be approximately inertially free above 
the pendulumn frequency. Nevertheless, thermal and other excitations of the suspension 
wires move the mirrors and lead to (narrow-band) data stream artifacts, located at the violin- 
mode resonance frequencies of the suspension wires. These narrow-band features, which 
stand some 40-60 dB above the shot noise background in the LIGO I detector, are signals 
themselves that, unless properly treated, can confound efforts to detect actual gravitational 
wave signals. 

In the waveband of the detectors now under development the transfer function, which 
determines how the detector output at a particular frequency is related to incident radiation 
at that frequency, is constant. The time evolution of the signal arising from any incident 
gravitational waves is thus determined solely by the waves themselves. The time evolution 
of the signal arising from suspension violin mode excitations, on the other hand, is strongly 
infiuenced by the the dynamics of the modes. We make use of that distinction to design a 
Kalman filter [5] that can, in real-time, detect violin mode excitations — thermal or other- 
wise — as distinct from actual gravitational wave signals. The measured excitation can be 
used to subtractively remove the modes from the detector output, significantly reducing the 
rms detector noise and simplifying the down-stream analysis. Additionally, these measured 
excitations can be used to diagnose the instrument's behavior. 

Other techniques have been described for the identification and removal of resonant line 
features in gravitational wave detector data. In particular, Sintes and Schutz [6,7] have 
discussed methods suitable for addressing nearly periodic artifacts that appear in a range 
of harmonics and whose phases, in each harmonic, are linearly related. These methods 
are not suitable for isolating the violin modes, however: though a full spectrum of violin 
mode harmonics do appear in the detector output, these modes are all independent and 
independently excited. Allen and Ottewill [8] have applied the multi-taper spectral method 
of Lees and Park [9] to estimate the contributions of these resonant features to the power 
spectrum. This non-parametric method does not and cannot provide a time-domain estimate 
of the mode state or its instantaneous contribution to the detector output and so cannot be 
used for conditioning the data prior to analysis. 

In section II we summarize the theory of Kalman Filtering, laying the ground for its 
application to identifying gravitational wave detector suspension violin modes. The con- 
struction of the Kalman filter for a particular problem requires a dynamical model for the 
underlying physical process giving rise to the signal. We discuss our model and how it comes 



to be represented as a set of state-space difference equations on the physical observables. 
In section V we investigate the performance of the fiher using both simulated data sets, 
whose "signal" and "noise" content is known exactly, and experimental data taken at the 
LIGO 40 M prototype detector during November 1994.^ We summarize our conclusions in 
section VI. 



II. THE KALMAN FILTER 
A. Introduction 

A Kalman filter is a mechanism for predicting the multi-dimensional state of a dynamical 
system from a multi-dimensional observable. The observable is assumed to be linearly related 
to the system's state, which is assumed to evolve linearly with time. The observable is 
assumed to be contaminated by an additive Gaussian measurement noise of known co- 
variance. The dynamical system described by the state may be driven by a known input in 
addition to Gaussian process noise of known co-variance. 

For example, the system may be a pendulum, the state may be the pendulum bob's 
instantaneous position and momentum, and the observable may be the bob's position. The 
process noise is the thermal fluctuation in the wire and the measurement noise in this case 
may be the photon shot noise background. 

When the process and measurement noise are Normal {i.e., Gaussian) the Kalman filter 
gives the optimal estimate of the system's state, in the least-squares sense. When the noise is 
more complex than Gaussian the Kalman filter gives the best estimate of the state assuming 
that only the second order statistical characterization of the noise is known. Alternatively, 
when the evolution equations are not linear or the observables are not related linearly to the 
state, then the Kalman filter has a natural generalization that makes it the optimal means 
of predicting the state from the observables in a least-squares sense [11]. 

The derivation of the Kalman Filter equations is straightforward and instructive.^ In 
the particular application dealt with in this paper we assume that there is no driving force 
(aside from the process noise) and that the system state and observables are needed only 
at discrete times; correspondingly, we limit our derivation to that case. Denote the state of 
the system at sample time k by il^lk] and the observable at sample time k by z[k]: 

I \u] — ( "^^^ state vector of the dynamical \ (o a \ 

y system at sample A;; a M x 1 vector J \ ■ ) 

\u\ — ( '^'^^ system observable at \ ro iv.'i 

^ ' \ sample k; a, N x 1 vector J ' 

The state i/? evolves from one sample to the next according to the linear equation 



^Preliminary results are presented in [10]. 

^For a more detailed discussion of Kalman Filters, see [5,12-14]. 



tl^lk + 1] = A ■ ^l}[k] + w[k], (2.1c) 

where A is a linear operator and w[k] is an instance of process noise: 

(Linear transformation that \ 
takes il^[k] to i^[k + 1]; (2. Id) 

a M X M matrix / 

(Process noise: a M-dimensional \ 
Normally distributed random (2-le) 

variable; a M x 1 vector. J 

The process noise is assumed to be zero mean and white. It is then fully characterized by 
its CO- variance matrix W: 

W := E[w ■ w'^], (2.1f) 

where E[-] represents an ensemble average across instances of its argument. The process 
noise co-variance is a symmetric M x M constant matrix. 

The observable z at sample k is related to the system state by the linear equation 

z[k] = C ■ il^[k] + v[k], (2.2) 

where G is a linear operator and v[k] is an instance of measurement noise: 

( The linear transformation \ 
C:= 



(2.3) 



that relates the state ■0[A;] 

to the observable z\k\] an 

\'N X M matrix ) 

(Measurement noise: a A^-dimensional 
Normally distributed random variable; | . (2-4) 

a A^ X 1 vector 

The measurement noise is assumed to be zero mean, white, ^ and have co- variance V: 

V:=E{v-v^\ (2.5) 

where i?[-] represents an ensemble average across instances of ^[A;] of the measurement noise. 
The measurement noise co-variance is a symmetric N x N matrix. 

Finally, we assume that the measurement noise and process noise are independent: 

E[v ■ w'^] = 0. (2.6) 



^If the measurement noise is not white, it can be whitened by a suitable, hnear redefinition of the 
observable. 



B. Predicting the state from observations 

Suppose that we have an estimate il)[k — 1] of the process state at sample k — 1. Denote 
the error in this estimate by e[k — 1]: 

e[k-l]:=il}[k-l]-^[k-l]. (2.7) 

This error is a random variable whose second moments we denote P[A; — 1]: 

P\k~l\:=E[e[k-l]-e^[k-l]]. (2.8) 

From the estimated state at sample A; — 1, we can form an a priori estimate of the state 
at sample k using the state transition operator A: 

^[k]:=A-^[k-l]. (2.9) 

Associated with il)[k] is an error e[A;] whose second moments we denote P[k\. 

e[k]:=il)[k]-^[k]. (2.10a) 

P[k]:=E\^[k]-e^[k]] (2.10b) 

= A-P[k-l\- A^ + W. (2.10c) 

The final equality follows from equations 2.10a, 2.1c and 2. If. 
From il)[k]^ we can make a prediction 5 of the observable z[k\. 

z[k] = C-^[k]. (2.11) 

Now compare the estimate 2 [A;] to the actual measurement z[k\. These will in general 
be different. The Kalman Filter uses the difference between the the actual and estimated 
measurement to generate a corrected estimate if) of the state: 

^[k] = i^[k] + K[k] ■ {z[k] - z[k]) . (2.12) 

The linear operator K[k] is the Kalman gain, which we calculate in section II C. The refined 
state estimate ify[k] has an associated error e[k] (cf. eq. 2.7) and covariance P[k] (cf. eq. 2.8): 

P={I-K C) P {IK Cf , (2.13) 

where / is the identity matrix and we have used 2.7, 2.10a and 2.12 to represent P in terms 
of the unknown K and the known C and P, all at sample k. 

C. Calculating the Kalman Gain 

Our object is to find the Kalman gain K that minimizes an appropriate measure of the 
error P in the state estimate ip at each sample k. 



Recall that P is a co- variance matrix of the error e; consequently, its eigenvectors describe 
linear combination of the state vector elements whose errors are uncorrelated. Correspond- 
ingly, the sum of the eigenvalues — the trace of P — measures the variance of the total 
error in the state estimate. Adopting tr P as a measure of the error in the state estimate, 
we choose the Kalman gain K to minimize trP: 

K[k] = P[k] -C^/ fv + C- P[k] ■ C^) , (2.14) 

where we have adopted the notation 

Q = S/R (2.15a) 

to denote Q as a solution to the linear equation 

S = Q R (2.15b) 

for known R, S. Correspondingly, equation 2.13 for P can be simplified: 

P = {I-K C) P {IK Cf 

= {I-K C) P-{I-K C) K {V + C K C^)K^ + K V K 

= {I-K C) P. (2.16) 

D. Summary 

To summarize, we begin with an estimate of the state t/) and the associated error P[k] 
at sample k — 1. From these, we 

1. Form an a priori estimate of the state ily[k], the error associated with that estimate 
P[k], the predicted measurement z[k], and the Kalman gain K[k], all at sample k (eqs. 
2.9, 2.10c, 2.11, and 2.14); 

2. Observe the system to determine z[k]; 

3. Form corrected estimates of the state iplk] and the estimated error P[k] (cf. eq. 2.12, 
2.16). 

We end as we begin, ready to make a new estimate of the state at sample k + 1. Once we 
are started, we can continue as long as observations can be made. To begin, we need only 
supply an initial estimate of the state and the error. These initial guesses are not critical as 
the algorithm continually corrects the estimated state and re-estimates the error with each 
iteration, converging on an optimum (in the sense of minimum tr P) estimate of the state. 
The cost per iteration — in either computation or memory — is constant. Only the 
present estimates of the state and the associated error need be retained from step to step: 
all of the history necessary to estimate the future values of the state and the associated error 
is encoded in these quantities. 
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III. APPLICATION: SUSPENSION VIOLIN MODES 
A. Introduction 

The cavity mirrors of terrestrial gravitational wave detector are suspended as pendula 
in order that they be inertially free at frequencies above the pendulum frequency. The 
suspension wires are under tension from the mirror mass and have a spectrum of "violin" 
modes. Each violin mode is, like the pendulum mode, endowed with ksT of energy from 
its contact with a thermal bath. The corresponding motion of the cavity mirrors leads to 
noise of thermal origin in the detector output. Other disturbances affect the suspension, 
increasing the overall noise power in the violin modes. 

To minimize the impact of the thermal noise on the detector's performance the suspen- 
sions are designed to have high Q. Laser shot noise dominates the detector noise power in a 
broad band about the violin mode resonant frequencies: only the peaks of the stochastically 
excited violin modes stand above this background. The peaks are strong, however: in the 
LIGO 40 M prototype they stand approximately 50 dB above the shot-noise floor and the 
contrast in the LIGO I detector will be approximately the same. In the LIGO 40 M proto- 
type the fundamental violin mode resonant frequencies are in the (570, 606) Hz band with 
Q's on order several times 10^, while in the LIGO I interferometer the fundamental violin 
modes will have frequencies near 330 Hz and Q's on order 10^ [15]. Figure 1 shows the power 
spectral density (PSD) of the differential-mode read-out (IF0_DMR0) of the LIGO 40 M pro- 
totype interferometer in the neighborhood of the violin mode resonances for a typical stretch 
of data taken during November 1994. The units of the PSD are ADC-counts^/Hz. The line 
at 600 Hz is the ninth harmonic of the omnipresent 60 Hz power-main. Note how the violin 
mode resonances stand-out as sharp peaks some 40 dB above an approximately white (over 
this bandwidth) background, as described. 

The combination of violin modes in an approximately white background maps naturally 
onto a model involving a stochastically driven dynamical process observed in the presence 
of white measurement noise: 

• Each violin mode corresponds to a distinct dynamical processes driven by an indepen- 
dent noise source (thermal or otherwise); 

• The observable is the motion of the cavity mirror as reflected in the gravity-wave 
channel of the interferometer; 

• The measurement noise is the laser shot-noise, and other technical noises that affect 
the gravity-wave channel but do not drive the violin modes. 

In the remainder of this section we describe how we restrict attention to the narrow band 
surrounding the violin mode, where the model associated with the use of a Kalman filter 
is particularly apt, and develop the state-space equations that describe the modes in that 
band. 



B. The Observations 

The violin mode resonances are very weakly damped; so, the noise power in the modes 
contributes significantly to the overall detector output only in a very narrow band about 
each mode's resonant frequency. Correspondingly, to monitor one or more violin modes we 
focus attention on a narrow band in the the neighborhood of the mode resonant frequencies. 
For definiteness, denote the gravitational wave channel at full bandwidth (sample rate fs) by 
g and let the band of interest, which includes the mode frequency /o, range from /c — A//2 
to fc + A//2, for A/ much less than fc- 

To exploit the narrow bandwidth of the violin modes we form the complex quantity 

g'[k] = exp (-27r^A;^') g[k]. (3.1) 

The spectrum of g' in the band [— A//2, A//2] is equal to the spectrum of g in the band 
[fc — A//2, fc + A//2]. Focus attention on just this narrow band by down-sampling g' to 
the new sample rate 2A/ using, e.g., a polyphase technique [16, Chapter 4].^ Denote the 
down-sampled g' by z, which is complex-valued. 

From z, it is convenient to construct the real vector sequence z[k]: 



^W= o.m • (3-2) 



^[z 

The vector- valued time series z[k] plays the role of the observation in our Kalman filter 
implementation. 

C. Dynamical Model 

The gravitational wave channel g is the sum of the violin modes (whose behavior we wish 
to determine) and other noise sources. In constructing z we focused attention on a narrow 
band about the mode resonant frequencies. Here we provide a model for the mode behavior 
as refiected in z. 

We first note that, near resonance, the details of the damping are unimportant; con- 
sequently, we can model each violin mode as an independent, viscously damped harmonic 
oscillator whose coordinate u is the amplitude of the mode present in the detector output. 
The coordinate u satisfies the differential equation 

u + ^u + u;^u = F{t), (3.3) 

where F is the (stochastic) driving force. The "gravity-wave" channel is the sum of u and 
other noises (and signals). 



^The down-sampling operation also involves the appropriate low-pass anti-alias filtering. 



We can relate u to the contribution of the viohn modes to z by mixing u with the phase 
factor exp{—iuct) (where Uc is equal to 2tt fc): 

ip = uexp{—iuct). (3.4) 

The complex ip satisfies the differential equation 

^0 , c^, . \ J. , / .2 2 , ,^c^0 



'^ + [~?T + 2ic^c j ^ + i^Q- ^c+ ^~7r- ) V' = Fexp{-iujct), (3.5a) 

or, in the Laplace domain, 

(s — p+) (s — p_) ■ijj{s) = F{s + iuc) (3.5b) 

where 



g{s) = / dt exp{—st) g{t) (3.5c) 



p^ = -^-^u.±^u.o^l-^. (3.5d) 

(Note that in equations 3.5 we have assumed that ip is causal — i.e., it vanishes at sufficiently 
earlier times — and used a two-sided Laplace transform. Written in this way the Laplace 
transform ip oi ip is related to its Fourier transform ip through the substitution of 2TTif for 
s: 

/oo 
dtexp{-2mft)git). (3.6) 

-oo 

(Note that we have adopted the engineering convention for the Fourier transform.) After 
down-sampling, the real and imaginary parts of ip are the contributions of the violin modes 
to the real and imaginary parts of z. 

The two poles p± in the Laplace transform solution ip{s) correspond to the positive and 
negative frequency resonances of the oscillator. When we choose Uc close to ujq we placed one 
of the resonances in ip near zero angular frequency and the other far away: correspondingly, 
one of the poles p± is near the origin and one is distant. Letting p_|_ be the "near" pole and 
P- be the "far" pole we write 

(.-pJV-W = ^<i±M. (3.7) 

s — p_ 

or, for \s\ <^ \p-\, 

{s-p+)iP{s) = -pZ^F{s + tuj,). (3.8) 

The spectrum of ip arising from this second equation, which differs from the first only by the 
source term, is identical to the spectrum of ip arising from the first as long as \uj\ is much 
less than |p_ | . Since down-samphng restricts the bandwidth of ip in just this way we adopt 
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equation 3.8 as the evolution equation for ip. Correspondingly, near zero angular frequency 
ip satisfies 

^ - p+^ = -pZ^F{t)e-''^^'. (3.9) 

Now write, instead of the complex ip, the real vector i/j whose components are the real 
and imaginary parts of ip: 

The first order differential equation 3.9 describing complex ip near zero frequency becomes 
a coupled set of first order differential equations 

'^ = A'-il^ + F (3.11a) 

where 



/ ^{pZ^e' 



■iuJct'\ 



^ = ^ ( aipiv- ) j ■ P.iic) 

Our observation is not of continuous z{t) but of sampled z[k]. Suppose that the final 
sampling rate is A/ and approximate the evolution of ij) as free between samples, with the 
accumulated effect of the stochastic force F acting at the sample times. Then the sampled 
if)[k] satisfies 

il)[k] = A ■ i/)[k - 1] + w[k] (3.12a) 

where 

A = exp(^^^ (3.12b) 

and w[k] is a random vector related to F. The quantity il^[k] in this case is equal to the 
contribution of the violin mode to the observable z[k], or 

z[k] = C ■ il^ik] + v[k], (3.12c) 

where v represents the (approximately) white background in which the violin modes are 
immersed. Equations 3.12 are the state-space form of the Kalman equations. To complete 
our description of the suspension modes in this form it remains to determine the properties 
of the process and measurement noises w and v. 

The statistical properties of w derive ultimately from the statistical properties of the 
force F (cf. 3.3). In applying the Kalman filter we assume that w, and consequently F, 
is white. Assuming that modes are driven principally by thermal forces we can use the 
Fluctuation-dissipation theorem [17] to find the power spectrum of the thermal force driving 
the mode: 
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Spif) = AksT^iZ) (3.13) 

where Z is the system's impedance. In our apphcation, 

Z = -^' + ^^--^0. (3.14) 

so, ^{Z) is independent of frequency, the thermal force acting on the modes is white, and 
the assumption made in applying the Kalman filter to this problem is satisfied exactly. 

The quantity v is the measurement noise: i.e., the noise within which the violin mode 
signal is embedded. In our application the measurement noise is associated with the laser 
shot noise. In the current generation of large detectors the shot noise has the spectrum 



Sh{f) = So 



1 + 



(3.15) 



with knee frequency fk approximately 100 Hz. the LIGO 40 M prototype the mode fun- 
damental frequencies are confined to an approximately 40 Hz bandwidth about 590 Hz; in 
the LIGO I detector they will be in a similarly narrow band about 350 Hz. As long as we 
are able to confine attention to bandwidths small compared to the mode frequency the shot 
noise in that band will be very nearly white; correspondingly, v will be white and have a 
covariance V proportional to the unity matrix. 

To summarize, we construct the observable z from the gravitational wave channel g 
by mixing it with a local oscillator and down-sampling to the bandwidth A/ about zero 
frequency. We model the sampled observable as arising from the equations 

tl^lk] = A ■ ip[k - I] + v[k - I] (3.16a) 

z[k] = C ■ il^[k] + w[k] (3.16b) 



where 



exp ( ^^ (3.16c) 



D. Reconstruction 

In applying the Kalman filter we determine, at each sample, the state estimate t/?[/c] from 
which we form an estimate z^ of the contribution of the violin mode to the detector output: 

z„ = C-'0. (3.17) 
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The violin mode contribution to the observation is an instrumental artifact that is unrelated 
to the incidence of gravitational waves on the detector; correspondingly, we may wish to use 
the prediction to subtractively remove the artifact from the data stream.^ 

The prediction z^ for the violin mode contribution is formed at reduced bandwidth 
and offset frequency. From z^ we can construct the predicted contribution violin mode 
contribution to the full bandwidth observation {g[k], cf. 3.1) by reversing the steps taken to 
form z from g (cf. sec. IllB). In particular, we first form the complex scalar sequence 

%={l i)-z, (3.18) 

and up-sample it to the full sample rate, forming G'^. The mixing operation undertaken in 
equation 3.1 can then be reversed, 

G^[k] = G'[k] exp (2mk^j . (3.19) 

The result is still complex: it is missing the contributions from negative frequencies that 
were discarded in the low-pass filtering that accompanied the original down-sampling oper- 
ation. Since the original observation was real, however, this contribution is just the complex 
conjugate of the positive frequency contribution, which we have; consequently, the estimated 
contribution of the violin mode to the data stream g[k] is 

g^k] = 2^{G[k]). (3.20) 

The residual difference 

gr[k]=g[k]-g4k] (3.21) 

is then free of violin mode artifact, to the accuracy of our estimation procedure. 

E. Estimating the measurement and process noise 

To use the Kalman filter we must have a model for the process dynamics (A) and how the 
process contributes to the observable (C). These we described for our particular application 
— the violin modes of an interferometer mirror suspensions (each mode an independent 
process) — in section IIIB. Additionally, we must have estimates for the process and 
measurement noise co- variance matrices {W and V). Here we describe how we estimate 
these from a short segment of filter input, taken at the beginning of the input. 

Consider the typical bandwidth of the observable z, shown in figure 1. The process 
contributes significantly over only a small fraction of the bandwidth, where the resonance 
rises above the measurement noise floor. To estimate the measurement noise, we form the 



^In section IV we show that the Kalman filter estimate of the violin mode contribution to the 
observation is insensitive to the presence of a gravitational wave signal. 
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PSD of z and take the average power far from the resonant peaks. This average power P^ 
is the total measurement noise (tr V). We assume that there is no cross-covariance between 
the noise in 3fJ(z) and '^i^z\ in which case 

V = ip„.(j;). (3,22) 

To estimate the process noise we focus on that part of the spectrum near the resonant 
peak where the process dominates. We assume, in our model, that the process is driven by 
white noise. From the Kalman filter model, the power spectral density of the contribution 
to z in this bandwidth is given by (cf. equation 3.9) 

P,{uj) = 2- (3-23) 

ui+(2Qu;,-uo{4Q^-lf') 

With the exception of the total process noise (trVF) the right hand side of equation 3.23 
is known. The left hand side we measure. To estimate tiW, we measure Pz and sum it 
over the limited bandwidth dominated by the process. Setting that equal to the same sum 
executed on the right hand side of equation 3.23 we determine the proportionality constant 
tr W. Assuming that the process noise is equally distributed between the real and imaginary 
parts of z, with no cross- variance, we thus determine the estimated process noise covariance: 

W=^-{tTW)(liy (3.24) 



IV. CHARACTERIZATION 

Before providing results from the application of our filter on LIGO 40 M prototype data 
we characterize the filter's behavior on "synthetic" or "mock" time series. These time series 
are constructed so that the violin mode state, process noise, measurement noise, and signal 
content are known exactly. By applying the filter to these data sets and investigating the 
difference between its estimation of these quantities and their actual, known values, we 
characterize the filter's performance. 

Four applications to synthetic data sets are described here: 

1. A single isolated violin mode, driven by Gaussian white process noise and observed in 
the presence of Gaussian white measurement noise (cf.sec. IV B); 

2. Several closely spaced violin modes, driven by non-Gaussian white process noise and 
observed in the presence of non-Gaussian white measurement noise (cf. sec. IV C); 

3. A single isolated violin mode, driven by Gaussian white process noise with a superposed 
impulsive excitation and observed in the presence of non-Gaussian white measurement 
noise (cf. sec. IV D); and 
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4. A single isolated violin mode, driven by non-Gaussian white process noise and observed 
in the presence of non-Gaussian white measurement noise with a superposed impulsive 
excitation (cf. sec. IV E). 

The first case demonstrates the filter operating in the regime where the measurement noise 
and process noise satisfy the assumptions made in the derivation of the Kalman filter. The 
second case examines the filter's ability to distinguish a violin mode from measurement noise, 
and to distinguish between two different modes, when the noise is not as simple as Gaussian 
and modes are closely spaced in frequency. In the third case we test the filter's ability to 
respond to impulsive excitations of the suspension modes. In the final case we introduce an 
impulse into the measurement noise and observe the filter's ability to discriminate between 
it and the contribution from the violin mode. 

As they appear in the detector data stream, violin modes are characterized entirely by 
their resonant frequency /o and damping constant Q. Since we are dealing always with 
discrete-time signal processing, the scale of the frequency does not matter: the only relevant 
quantity (aside from Q) is the ratio of the resonant frequency to the sampling frequency fs- 
Nevertheless, to foreshadow our application of the filter to data taken at the LIGO 40 M 
prototype, we assume that the sampling frequency is 9868.421 Hz, the sampling rate of the 
LIGO 40 M prototype, and we always choose the simulated violin modes to have resonant 
frequencies coincident with the frequencies of modes identified in that instrument. With 
this convention, the frequency and damping constants of the simulated modes used in the 
examples to follow are given in table I. 

A. Generating synthetic data sets 

Before discussing the results of our investigations into the Kalman filter's behavior we 
describe how we generated the random time series used in those investigations. 

In all cases the measurement noise is white: i.e., its power spectrum was constant over 
its bandwidth. For Gaussian noise measurement noise we used Matlab's [18] randn random 
number generator. Non-Gaussian noise was modeled by a two-component mixture- Gaussian 
model: 

P{x\a, Hi, H2, cTi, 0-2) = aN{Hi, ai) + (1 - a) A^(/i2, cr2) (4.1a) 



(4.1b) 



Here N{fi, a) represents a Normal process with mean fi and variance a^, the parameters 
a, Hk and a^ are fixed with < a < 1. Figure 2 shows the distribution of this mixture 
Gaussian process together with the distribution of a Gaussian process with the same mean 
and variance. The mixture Gaussian distribution described by the parameters in equation 
4.1 is strongly leptokurtic [i.e., it has a significant excess of high amplitude events). 

To generate the contribution to the observation from the process we begin with the 
differential equation modeling the process (cf. 3.3): 
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a = 0.7, 






Hi = 0, 


H2 = 


= 0, 


ai = 1, 


0-2 = 


= 1.5. 



r + ^f + ujlr = F, (4.2) 

where r is the position of the suspended mass. We make no assumptions here about driving 
force F. The detector gravitational wave channel consists of a superposition of measurement 
noise v and the oscillator coordinate r. The solution to this equation can be expressed as 
an all-pole linear filter H acting on the driving force F. In the Laplace domain the filter is 
given by 

f{s) = H{s)F{s) (4.3a) 



where 



s^ + suoq/Q + uol 
(s-p+)(s-p_) 



p± = - {u;o/2Q ± luo Vl - l/4g2J (4.3c) 



The frequency response of the filter at angular frequency u is equal to H{s) evaluated on 
the imaginary axis: 

H{lu) = H{iuo). (4.4) 

The two-pole filter H operates in continuous time while the detector output we wish to 
simulate is discretely sampled at a rate /s, which we assume is at least twice as great as the 
resonant frequency /q. Using the analog filter H{s) as a prototype, we design a digital filter 
with a similar response at frequencies much less than the Nyquist frequency. 

Suppose that the discretely sampled sequence r[A;], driven by the input F[k], satisfies 

4k] = Y,mnk-J] -J2a[j]u[k-j], (4.5) 

j=0 j=l 

where a[j] and b[j] are the constant filter coefficients. Define the z-transform of the sequence 

y[k] by 

oo 
fc=— oo 

(The z-transform of a digital filter is analogous to the Laplace transform of a continuous 
filter.) In 2;-space, the hnear filter described in equation 4.5 can be written 



^A^, 



EZom^- 



■3 



f{z) = ^'=1 '■"' F{z) (4.7a) 

G{z)F{z) (4.7b) 
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The frequency response of the digital filter G is equal to G{z) evaluated on the unit circle 
in the complex z-plane: 

G {e^^'f/f') (4.8) 

where /^ is the sample rate of the filtered sequence. If, by a conformal transformation, we 
map the right half s-plane on to the interior of the unit circle in the z-plane, then we will 
have mapped the analog prototype H{s) on to the digital filter G{z): 

G{Z) = B{s)\s=2Q{z-l)/(z+l), (4.9) 

where Q is the so-called warping constant. This mapping will not be without distortion; 
nevertheless, at frequencies much less than the Nyquist frequency the distortion will be 
minimal. We can choose the warping constant Q to guarantee that the digital filter's response 
matches exactly the analog filter's response at exactly one frequency: 

G{e-^^'f'/f^) = H{2nif') (4.10) 

when 

71 f 

= . (4.11) 

tan(7r/V/,) ^ > 

We choose the warping constant so that the response of the digital filter matches the response 
of the analog filter at the resonant frequency /q. This completes the design of the digital 
filter G{z). Given a pseudo-random process F[k] and the recursion relationship of equation 
4.5 we evaluate r[k], the simulated contribution of the violin mode to the detector output. 
Finally, to obtain the simulated gravitational- wave channel output, we add to r[k] another 
pseudo-random process n[k], which represents the measurement noise. 

In our simulations we used the esc Matlah [18] Normal random number generatorrandn 
wherever Normal deviates were called for: i.e., in forming either the measurement noise, the 
process noise or the stochastic driving force F . 



B. Single mode: Normal, w^hite process and measurement noise 

Using the techniques described above we generated several hundred seconds of simulated 
data, consisting of a single mode with the frequency and damping constants given in the 
first column of table I. The process noise was white with variance of 2.38; the measurement 
noise was also white with variance 2.8. 

After generating the data set we applied the Kalman filter, focusing on a 6 Hz bandwidth 
A/ centered about the mode's resonant frequency. Following the discussion in section III E 
we estimated W and V using approximately 100 s of the simulated data.^ We applied the 



^Since this data set is simulated the values of W and V are known a priori; however, in a real 
application they will not be known. We used the same procedure to estimate the filter as would 
be used in the field. 
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Kalman filter to the subsequent simulated data, using these values of W, V, setting the 
initial state estimate i/' to zero and the initial estimated error P equal to WI (where I is 
the identity matrix). Figure 3 summarizes, in three panels, the Kalman filter's effectiveness 
in estimating the mode state. The upper panel shows the power spectrum of the simulated 
observations in a wide band about the mode frequency. The middle panel shows the power 
spectrum of the Kalman estimate of the mode contribution to the observable. Finally, the 
lower panel shows the power spectral density of the difference between the actual data and 
the estimated contribution from the mode.^. Comparing these three panels it is clear that 
the Kalman filter has correctly identified the modes contribution to the observable. 

Since the data is simulated we know exactly the measurement noise and mode contri- 
butions to the observable. Correspondingly, we can evaluate the actual error in the filter's 
estimate of the mode state. The mode (estimated) state is described by the complex ('^/'[A;]) 
ip[k], which we describe by its amplitude and phase 

A := \iIj\ (4.12a) 

tan0:=S(V^)/3f?(^). (4.12b) 

The difference between the actual and estimated mode states we describe by the quantity 

A: 



A := ^-!-^ — Z^^. (4.13) 

\m+m\ 

Figure 4 shows the instantaneous estimated mode amplitude, phase and error A. The top 
panel shows the estimated mode amplitude. The middle panel shows the estimate mode 
phase, less 27r/ot, where /q is the mode frequency. The bottom panel shows the error in the 
estimated mode state. Following an initial epoch, during which the filter is "locking-on" to 
the mode state, the relative error in the state estimate falls quickly to zero. 

The middle and bottom panel of figure 3 show a slight depression in the power spectral 
density of the residual in the band that the Kalman filter operates. This slight depression 
is associated with the error in the state estimate, which is shown in the bottom panel in 
figure 4 to be on order several percent. The depression is uniform throughout the band, 
showing that the errors introduced by the subtraction of the estimated mode contribution 
from the data are white. For this reason it is unlikely that these errors can be reduced by, 
for example, a "better" estimate of W, V, /o or Q: changing any of these would lead to an 
error that was not white, but reflected either an over- or under-estimate of the power in the 
mode relative to the measurement error. 

We expect that the continuously updated mode amplitude and phase estimates as a 
function of time, either analyzed directly or displayed in a strip chart of which the top two 
panels of flgure 4 may be considered a snapshot, will be particularly useful for monitoring 
the state of the detector. 



^While the power spectral density in a relatively wide band about the mode resonant frequency 
is shown, subtraction of the estimate affects only a 6 Hz bandwidth about the mixing frequency /c 
(cf. sec. IIIB and eqns. 3.4, 3.19). 
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C. Multiple modes; non-Gaussian white measurement and process noise 

In deriving the Kalman filter we assumed that the process and measurement noise were 
white; however, we never required that they be Gaussian. Here we investigate the effec- 
tiveness of the filter when these noises are non-Gaussian and when there are several modes 
present, at least two of which are very close together in frequency. This is a critical appli- 
cation, since in actual data there are multiple closely spaced violin modes resulting from 
several suspension wires of the mirrors, which cannot be estimated separately (cf. 1). The 
mode frequencies and damping constants for this example are given in the second column 
of table I. 

For isolated modes {i.e., those whose resonant peaks do not overlap with each other 
above the level of the measurement noise), one can design and apply separately a Kalman 
filter for each mode. Closely spaced modes cannot be treated in isolation; however, since the 
mode dynamics are independent, the development of a Kalman filter that treats N modes 
jointly is straightforward: 

• Referring to section III B, fix a carrier frequency fc and a bandwidth A/ such that the 
modes of interest all have resonances in the band [fc — A//2, fc + A//2]. 

• Considered separately, the dynamics of the mode J (1 < J < A^) and the contribution 
Zj of the mode to the measurement are described by the equations 



z,[k] = Cr^,[k]. 



(4.14a) 
(4.14b) 



The Kalman equations describing the joint evolution of the A^ modes, and their con- 
tribution to the observation z are then 



^[k] = A-^[k-l] +w[k-l] 
z[k] = C ■ ^[k] + v[k], 



(4.15a) 
(4.15b) 



where 



^ := 



f^i\ 



S(^i 



3fJ(^7v) 
(A, 



\ 



An J 



(4.15c) 



(4.15d) 



(4.15e) 





(Wi 




W: = 








\ 


Wn 




(V, 


\ 


V: = 








\ 


Vn / 


c--={ 


Ci ■■ 


■ Cn) 




^1 ■ 


■■ 1 0\ 




,01- 


■■ 1 j 



(4.15f) 

(4.15g) 

(4.15h) 
(4.15i) 

and W j and V j are, as before, the process and measurement noise. 

In constructing W we have approximated the individual modes as being driven by inde- 
pendent process noise sources. Correlated process noise leads to correlations in the states of 
the individual modes. Even in the absence of an explicit correlation in the assumed process 
noise covariance {i.e., off-diagonal terms in W), the filter will, if warranted by the observa- 
tions, correlate the state estimates. Lacking knowledge of the driving noise (beyond contact 
with independent thermal baths), assuming that the process noise is uncorrelated guaran- 
tees that there is no prejudice that the mode states are correlated without constraining the 
state estimates to be uncorrelated. 

To evaluate the effectiveness of the Kalman filter in treating multiple modes driven by 
non-Normal noise, we applied a single Kalman filter to this three mode system, choosing 
/c equal to 578.6 Hz and A/ equal to 18 Hz. The top panel of figure 5 shows the power 
spectral density of the simulated gravity wave channel in a 60 Hz bandwidth about the 
carrier frequency /c, while the bottom panel shows the power spectral density of the residual 
difference between the simulated observation and the filter's prediction of the contribution 
from the violin modes. The residual broad feature that remains following the subtraction 
of the mode estimate is due to the non-Gaussian nature of the process noise that drives the 
mode. The filter has made a best-estimate of the mode state, given only knowledge of the 
driving force variance. In this example, however, the variance provides a very incomplete 
picture of the process noise. Still, the filter has clearly identified approximately 45 dB of 
the approximately 55 dB of the mode contribution to the measurement. 

The state ^ is the direct sum of the state vector estimates i/'j of the individual modes; 
consequently, the filter identifies not just the total contribution of each mode to the observa- 
tion, but the state of each mode separately. The four panels of figure 6 show the estimated 
amplitude and phase (less 27r/jt for each mode) of each of the two overlapping modes at 
584.6 Hz (/2) and 585.4 Hz (/s). In this example the two modes are known to be uncorre- 
lated; consequently, we expect that the state estimates should also be uncorrelated, which 
we have found to be the case. 
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D. Effect of a transient excitation of a mode 

To evaluate the effectiveness of the Kalman fiher in estimating the filter state in the 
presence of a transient excitation, we have simulated a single mode observed in the presence 
of the same mixture-Gaussian measurement noise described above (cf. equation 4.1). The 
process noise is taken to be Gaussian with variance 2.38. At 20 seconds into the simulation 
we add to the process noise a single pulse of magnitude 35o", where a is the rms amplitude 
of the process noise. The top panel of figure 7 shows the estimated amplitude of the mode 
while the bottom panel shows error A in the mode state estimate: 

A := ^I^W-fWI (4.16) 

\i'\k]+i\k]\ 

After the filter's start-up transient the error in the estimated state is small until the transient 
excitation. There is, not surprisingly, a sudden increase in the error at the moment of the 
excitation; however, the filter quickly adjusts to the new state. 



E. Effect of transients in the measurement noise 

Sudden excitations in the observed time series may arise from sources other than excita- 
tions of the violin modes. To evaluate the effectiveness of the Kalman filter in discriminating 
between excitations of the mode and measurement noise transients we simulate a single line, 
driven by Gaussian noise. This mode is embedded in the same mixture-Gaussian measure- 
ment noise as described in section IV C, but with a pulse in the measurement noise 14 
seconds into the simulation. The pulse lasts for a single sample and has an amplitude of 
350", where a rms amplitude of the measurement noise. 

The lower panel of figure 8 shows A (cf. 4.16), measuring the magnitude of the difference 
in the estimated and actual mode state; the upper panel shows the estimated state itself. 
The error in the mode state estimate always remains quite small, showing that the filter 
discriminates between excitations of the modes and excitations in the measurement noise. 

F. Summary 

The Kalman filter described here is able to correctly determine instantaneous state of 
the violin mode suspensions. This is true whether the mode is driven by, or observed 
in the presence of, Gaussian or non-Gaussian noise. The filter does not require that the 
measurement or process noise be stationary: it correctly tracks the mode after it has been 
impulsively excited and can distinguish between excitations of the mode and transients in 
the measurement noise. 



V. VIOLIN MODES IN LIGO 40 M PROTOTYPE DATA 

In this section we discuss the application of the Kalman filter to data taken at the 
LIGO 40 M prototype detector during November 1994. The general parameters of this run 
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are described elsewhere [19]. The particular results presented here all refer to the particular 
locked segment beginning at GPS time 469357019 s on 19 November 1994 and lasting for 
2665.54496 s. 



A. General character of the resonant features 

At the time of the November 1994 data run each of the four cavity mirrors in the 
40 M prototype was suspended by four cylindrical wires. Setting aside the coupling of the 
wires through their attachment to the suspended mirror there are two fundamental modes 
associated with each wire: one with motion nominally along the optical axis and one with 
motion perpendicular to that axis. Owing to imperfections in the wires, the attachments 
of the wires to the test masses, and the coupling of the modes through the motion of the 
mass all the modes will have some component along the optical axis. We thus expect 32 
fundamental violin modes. 

A high resolution power spectrum focused on the bandwidth (570 Hz, 606 Hz) and taken 
of the data from a single locked-segment {i.e., a single epoch during which the interferometer 
operated without interruption as a gravitational wave detector) of 44 m duration shows 28 
distinct resonant features. One of these is clearly identified as the ninth harmonic of the 
60 Hz power-main, leaving 27 different modes whose frequencies and quality factors are 
listed in table HI. 

Referring to table HI note that, apart from modes 7, 16, 19, 24, and 25, the lines appear 
in pairs, with a high amplitude line closely associated with a lower amplitude line at a 
slightly higher frequency (less than or of order 0.2 Hz). This association of lines is clearly 
shown in figure 9. Of the five apparently unpaired lines, three show companions in a high 
resolution power spectrum taken in the band of the first harmonic, suggesting that all five 
of the "missing modes" may be so closely spaced in frequency with other, identified modes, 
that we are unable to resolve them. 



B. Apphcation 

Using the techniques described in section HI we have applied a Kalman filter designed to 
track the 27 violin modes identified in table HI. The top panel of figure 10 shows the power 
spectrum of the IFO_DMRO (interferometer differential mode read-out, or "gravity wave") 
channel in the 45 Hz band beginning at 565.0 Hz, with the violin modes clearly present. 
The bottom panel shows the power spectrum of the residual after the modes identified by 
the filter have been subtractively removed. Very little power associated with these modes 
remains.^ 

Figure 10 also shows an overall approximately 2 dB reduction in the background noise 
level within the band, which is due to the "wings" of the predicted contribution of the violin 
modes to the total noise. Our experience with simulated data (cf. sec. IVB-IVE) shows 



^The 600 Hz line, which is the ninth harmonic of the power-main, was not modeled in our analysis. 
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that the filter generally makes fractional prediction errors whose magnitude is on order 1% 
and that, as we have implemented the filter, this error tends to systematically over-estimate 
the contribution of the mode wings to the total noise. This is consistent with spectrum of 
the residual found here. 



C. Preliminary characterization of mode motion and residual 

The Kalman filter provides us with a decomposition of the observation into the contri- 
butions from the individual violin modes and a residual, which is, in a least-squares sense, 
as free from the influence of the modes as we can make it. The noise character of the vio- 
lin modes, which are instrumental artifacts, can thus be studied separately from the noise 
character of the residual. 

Figure 11 shows in histogram form the frequency with which samples of given ampli- 
tude (normalized by rms amplitude) appear in the (570, 595) Hz bandwidth of the original 
data ("raw"), in the contribution to this band from the modes as identified by the Kalman 
filter, and in the residual difference between the raw data stream and the estimated con- 
tribution from the modes. If the detector noise were strictly Gaussian the distribution of 
squared- amplitudes (normalized by the mean-squared amplitude) would follow an exponen- 
tial distribution and the resulting curves would be, on the semi-log graph as shown, a straight 
lines with slope — log^g e. We have evaluated x^ for each histogram, using the left-most 47 
bins shown in figure 11 and with the remainder of the data bins collapsed into a single 
additional bin, so that no bin has fewer than 40 events. The same binning was used for 
all three x^ statistic evaluations. Table IV gives the resulting x^, the number of degrees of 
freedom, and the corresponding p-value. 

Referring to table IV we see that the Kalman filter has resolved the total noise in this 
band into two components that have very different noise character. While in all cases the 
statistical character of the noise is far from Gaussian, the contribution identified as arising 
from the violin modes is much better behaved {i.e., much closer to Gaussian) than the 
residual. The close correspondence between the x^ for the "raw" data and the identified 
violin mode contribution reflects the fact that the violin mode contribution, even while it is 
narrow-band, is of such great amplitude that it dominates the total detector noise in this 
band. Recalling our earlier results (cf. sec. IV D and IV E), which show how the Kalman 
filter is able to discriminate correctly between excitations of the modes and impulses in the 
measurement noise, we conclude that the difference in the noise character is real and that, 
in the LIGO 40 M, the violin mode contribution to the overall noise has many fewer outliers 
than the non-mode background. Since, in a real detector, it is the residual that will contain 
the gravitational wave signal, the importance of removing instrumental artifacts like the 
violin modes in order to discover the noise character of the residual cannot be overstated. 

A still closer examination, focusing on the contributions from the individual modes, 
shows that certain modes have much greater excess noise than others. Figure 12-14 shows, 
for twelve different violin modes, the frequency distribution of mode amplitudes. All of the 
modes shown here have a strong excess of high amplitude noise events: i.e., they are subject 
to strong, non-thermal forces. 
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VI. CONCLUSIONS 

Violin modes of the test mass suspensions are a significant instrumental artifact in the 
signal band of modern interferometric gravitational wave detectors: at the mode resonant 
frequencies they stand between 40 and 50 dB above the ambient noise floor. These modes 
are not excited or otherwise disturbed by the passage of a gravitational wave; yet, they 
contribute almost all of the noise power in a wide bandwidth. Correspondingly, the ability 
to identify and remove their contribution to the detector noise and characterize the residual 
noise can significantly increase the detector's sensitivity in the affected wave-band. 

Treating the modes as a stochastic signal in a measurement noise background we have 
developed and demonstrated the use of a Kalman filter for identifying the state of the 
suspension wires and the contribution of the violin modes to the detector noise. A Kalman 
filter uses the known dynamics of the process that underlies a noise component (in this 
case the stochastic excitation of the suspension violin modes) to identify that component's 
contribution to total system noise from measurements made on the total noise alone. When 
the stochastic driving force is Normal {i.e., Gaussian) the Kalman filter is optimal in the 
least squares sense. When the process noise is non-Gaussian, the Kalman filter estimate 
is optimal when only the second order statistics of the noise {i.e., its power spectrum) is 
known. 

Using simulated data, whose mode and measurement noise contribution is known exactly, 
we have shown that the filter works well in the presence of both Gaussian and non-Gaussian 
measurement noise and excitations of the violin modes. We have also shown that the filter 
easily discriminates between transients in the measurement noise and transient excitations 
of the violin modes themselves, capturing the latter in its estimate of the mode state while 
rejecting the former as not arising from excitation of the mode. Using data taken at the 
LIGO 40 M we have demonstrated the use of the filter by isolating and subtractively remov- 
ing the noise owing to 27 distinct violin modes, each standing between 40 and 50 dB above 
the ambient noise floor, over an approximately 40 Hz bandwidth. 

The Kalman filter describes here operates in the time domain, providing instantaneous 
estimates of the suspension violin modes states. The computational cost is negligible, making 
it suitable for use as an on-line diagnostic as well as for on- or off-line data analysis. 

Thermal and other excitations of the violin modes dominate the detector noise over 
a wide band. This additive contribution to the detector output will never include any 
gravitational wave signal, however. Separating the mode contribution from the total noise 
dramatically increases the detector's sensitivity by significantly reducing the mean-square 
noise amplitude. It also, however, reveals a new "layer" of non-Gaussian and transient noise 
that must be characterized and, where possible, removed through its association with other 
instrumental and physical environment monitors. Using the Kalman filter to separate the 
total noise into a violin mode contribution and a residual we find that, in the LIGO 40 M 
prototype, the noise owing to the violin modes is much "cleaner" {i.e., much closer to 
Normal) than the residual noise, which is strongly leptokurtic. Since, in a real detector, it 
is the residual that will contain the gravitational wave signal, the importance of removing 
instrumental artifacts like the violin modes in order to discover the noise character of the 
residual cannot be overstated. As we dig deeper into the detector noise the challenges of 
identifying the origin of the newly uncovered noise components becomes greater; however, 
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this is a challenge that must be faced and overcome if we are to make the most of the 
opportunities provided by the these new detectors. 
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TABLES 



Case 1 



Case 2 



Case 3 



Case 4 



/o (Qo) 

h [Qi) 

/2 (Q2) 



571.6 Hz (57 

N/A 

N/A 



10^) 571.6 Hz (57 x 10^) 
584.6 Hz(57 x 10^) 
585.4 Hz(57 x 10^) 



571.6 Hz (57 

N/A 
N/A 



10^) 571.6 Hz (57 x 10^) 

N/A 
N/A 



TABLE L Constants characterizing the simulations described in section IIIB. The sample rate 
in all cases was 9868.421 Hz. 



Line 


a 


/^i 


c^i 


^2 


^2 


(Hz) 












571.6 


0.7 





1 





1.5 


585.0 


0.7 





1 





1.5 


585.6 


0.7 





1 





1.5 



TABLE n. The mixture-Gaussian model parameters used to describe the non-Gaussian process 
and measurement noise used in section IV C. 
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Line 


Frequency 


FWHM 


ID# 


(Hz) 


(Hz/100) 


1 


571.5869 


1.26 


2 


571.6820 


2.10 


3 


578.3400 


1.01 


4 


578.4185 


1.40 


5 


578.6820 


2.00 


6 


578.8050 


1.31 


7 


581.0550 


2.00 


8 


582.3957 


1.70 


9 


582.5586 


2.60 


10 


583.5750 


1.55 


11 


583.7379 


1.64 


12 


583.9429 


1.90 


13 


584.1077 


2.20 


14 


594.1913 


1.70 


15 


594.2902 


2.20 


16 


595.2734 


1.30 


17 


595.9235 


1.10 


18 


596.1029 


1.50 


19 


597.6456 


1.80 


20 


597.9212 


1.50 


21 


598.1042 


1.40 


22 


598.9301 


2.00 


23 


599.0271 


1.30 


24 


599.1425 


1.40 


25 


599.3713 


1.18 


26 


605.3880 


4.50 


27 


605.5310 


1.20 



Q Amp. 

(10^) m_ 

4.55 58.0 

2.72 38.0 

5.73 68.0 
4.13 43.0 
2.89 56.0 
4.42 46.0 

4.42 57.0 

3.43 58.0 
2.24 57.0 
3.76 60.0 

3.56 40.0 

3.07 57.0 
2.66 38.0 
3.50 58.0 
2.70 42.5 
4.58 61.0 
5.42 60.0 
3.97 46.0 
3.32 58.5 
3.99 57.0 

4.60 40.0 
2.99 60.0 

4.61 36.0 
4.28 60.0 

5.08 47.5 
1.34 55.0 
5.05 45.0 



TABLE HI. Violin mode frequencies, full- width at half maximum (FWHM), derived damping 
constants Q, and mode amplitudes for the 27 violin modes found observed IFO_DMRO data channel 
of the LIGO 40 M prototype. The modes and their characteristics were determined from high 
resolution power spectra formed from the data sets themselves. There was no significant variation 
in the mode frequency or quality over the three days covered by the November 1994 data set. 
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total violin residual 

noise mode 



X^ 39.8514 39.7136 47.1459 

iVfrccdom 47 47 47 

p-value 0.7607 0.7656 0.4666 

TABLE IV. The correspondence between the statistical character of the LIGO 40 M detector 
noise, in the band 570-595 Hz about the violin modes, and a Normal distribution. The figure of 
merit used in the comparison is x? ■ The first column gives x^ for the total noise, the second column 
for that part of the noise attributed to the violin mode and the third column for the difference 
between the total noise and the violin mode contribution. The violin mode contribution dominates 
the total noise in this band; when it is removed, the residual is seen to be very poorly behaved. 
For more detail see the discussion in section VC. 
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FIGURES 



565 



Caltech 40m data, 19Nov94.4, Locked004.nc,44.426 minutes, f =9868.421 Hz 
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FIG. 1. Typical power spectrum of the IFO_DMRO channel of the LIGO 40 M prototype, showing 

the violin modes between 571 and 606 Hz. The units of the power spectrum are (ADC counts)^ 

Hz"^. Some lines are so closely spaced that they are not resolved on this figure. The line at 600 Hz 

is the ninth harmonic of the 60 Hz power-main. 
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solid:mixture Gaussian a=.7,o^=1,a2=1.5; dashed: Gaussian witli same mean and variance 




-4-3-2-1 1 2 3 4 

Amplitude 
FIG. 2. To test the effectiveness of the Kalinan filter when the measurement and process 

noise is not Gaussian we used simulated data drawn from a strongly leptokurtic mixture Gaussian 

distribution. The solid curve shows the non-Gaussian noise distribution while the dashed curve 

shows, for comparison, a Gaussian distribution of the same mean and variance. 
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FIG. 3. The effectiveness of the Kalman filter in identifying the violin mode contribution to 

the signal can be seen by comparing the power spectrum of the observed time series (top panel), 

the power spectrum of the Kalman estimated contribution of the mode to the observation (middle 

panel) and the power spectrum of observed time series less the filter's prediction of the violin 

mode's contribution to the observation (bottom panel). The bottom panel shows that the residual 

is, as expected, very nearly white. Kalman filter is applied to only a 6 Hz bandwidth about the 

line; this is the reason for the rapid fall-off of the predicted contribution of the mode to the total 

noise. The slight depression in the residual noise level in the signal band is discussed, along with 

other details, in section IV B. 
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100 120 

Time (s) 
FIG. 4. The Kalman filter produces instantaneous estimates of the mode state. Here we rep- 
resent those estimates as an ampUtude (top panel) and phase (relative to the resonant frequency 
angular frequency 27r/ot; center panel). In this example we know the actual state as a function of 
time. The bottom panel shows the instantaneous error in the state vector relative to its instanta- 
neous magnitude. For more details see section IV B. 
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FIG. 5. The effectiveness of the Kalman filter in identifying separately the contribution from 

two closely spaced violin modes can be seen by comparing the power spectrum of the observed 
time series (top panel) and the observed time series less the filter's prediction of the mode con- 
tributions. In this example the process noise and the measurement noise are both taken to be 
strongly leptokurtic. Despite the strong non-Gaussian character of the process and measurement 
noise the filter is able to identify all but 3 x 10~^ of the power in the mode. For more details see 
section IV C. 
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FIG. 6. The Kalnian estimates of the states of the two closely spaced lines in the example of 
section IV C. The right panels show the estimated amplitude (top), phase (less 27r/of; middle) and 
fractional error A (bottom) for the /o = 585.4 Hz line; the left panels show the same quantities 
for the /o = 584.6 Hz line. As is apparent from the errors the filter has successfully discriminated 
between the two lines. 
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FIG. 7. Can the Kalman filter track a mode that undergoes an impulsive excitation? At 20s 

into the simulation the mode was excited by a large force. The top panel shows the amplitude 

of the estimated mode state; the bottom panel shows the fractional error A (cf. eq. 4.16) in the 

mode state estimate. The estimate is momentarily upset; however, the filter quickly re-acquires 

the actual mode state. For more detail see section IV D. 
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FIG. 8. Can the Kalman filter discriminate between an impulsive excitation in the measurement 

noise and the contribution of the mode? A mode, stochastically driven by Gaussian process noise, 

is embedded in a mixture-Gaussian measurement noise. At 14s a large transient is added to the 

process noise. The top panel shows the Kalman estimate of the mode state amplitude throughout 

the simulation. The bottom panel shows the fractional error A (cf. eq. 4.16) in the estimated state. 

There is no evidence that the filter is sensitive to the transient measurement noise. For more detail 

see section IV E. 
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19nov94.4/Locked004.nc, 20.43 minutes, f =9868.421 Hz. 
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FIG. 9. A high resolution power spectrum of the data set shows that most of the large amplitude 

violin mode resonances have a lower amplitude satellite line at slightly higher frequency. Here we 

shows the high amplitude lines together with their satellites in the 583.3 and 584.3 Hz band. For 

more detail see section VA. 
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Caltech 40m data, 19nov94.4/Locked004.nc, 20.43 minutes, f =9868.421 Hz. 
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FIG. 10. Top panel: A high resolution power spectrum of the LIGO 40 M prototype data, 

focusing on the 571 to 605 Hz band, shows the violin modes standing 40 to 50 dB above the noise 

floor. Bottom panel: A power spectrum of the same data, after the Kalman filter prediction for 

the mode contribution to the measurement has been subtractively removed, shows that the modes 

have been removed, leaving only a small residual in addition to the measurement noise. For more 

details see section VB. 
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FIG. 11. Histograms of squared-amplitudes in units of the amplitude variance. The top panel 
shows the "raw" data from the IFO_DMRO channel in the 570-595 Hz band. The middle panel 
histograms the amplitude of the Kalman estimate of the violin mode contributions to IFO_DMRO, 
while the bottom panel histograms the amplitude of the difference between the raw data and the 
Kalman estimate of the mode contributions. In all cases the straight line is the expected result for 
Gaussian noise. The raw data "breaks" at an amplitude-squared S/N of approximately 6 (which 
corresponds to the level at which single events are expected in each bin), after which it shows 
clearly an excess of high amplitude events. The Kalman estimate of the contribution by the violin 
mode lines shows a very similar character, supporting the conclusion that the violin modes are 
dominating the noise in this band. The residual, on the other hand, appears to show a modest 
break at an amplitude-squared S/N of 4 and then a dramatically larger number of high amplitude 
events above a S/N of 6. For more discussion see section VC. 
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FIG. 12. Histograms of the square of Kalman estimated amplitude divided by mean amplitude 
for each of the individual Kalman estimated violin modes. For more discussion see section VC. 
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FIG. 13. Histograms of the square of Kalman estimated amplitude divided by mean amplitude 
for each of the individual Kalman estimated violin modes. For more discussion see section VC. 
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FIG. 14. Histograms of the square of Kalman estimated amplitude divided by mean amplitude 
for each of the individual Kalman estimated violin modes. For more discussion see section VC. 
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